Benard-Marangoni Convection in Two Layered Liquids 
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We describe experiments on Benard-Marangoni convection in horizontal layers of two immiscible 
liquids. Unlike previous experiments, which used gases as the upper fluid, we find a square planform 
close to onset which undergoes a secondary bifurcation to rolls at higher temperature differences. 
The scale of the convection pattern is that of the thinner lower fluid layer for which buoyancy and 
surface tension forces are comparable. The wavenumber of the pattern near onset agrees with the 
prediction of the linear stability analysis for the full two-layer problem. The square planform is 
in qualitative agreement with recent one- and two-layer nonlinear theories, which fail however to 
predict the transition to rolls. 
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Convection in fluids has been a fruitful system for the 
study of nonlinear, noncquilibrium patterns for almost 
100 years. QJ^H The original experiments of Benard 
|jj used shallow layers of whale oil, heated from below 
and open to the air above. Many years passed before it 
was conclusively shown that surface tension gradients, or 
"Marangoni", forces were crucial. In general, both 
Marangoni and buoyancy forces are present. HQ Begin- 
ning with Pearson || , theories of this instability have tra- 
ditionally neglected the dynamics of the upper fluid (air 
in Benard's experiments), replacing it by a constant heat 
flux boundary condition on the lower fluid. This con- 
dition is experimentally unrealizable, however. In well- 
controlled experiments p|,p[ pT|JTl| , p^ , p^[ | , the upper fluid 
is bounded above by a plate, on which a constant tem- 
perature boundary condition is maintained. The heat 
flux across the interface is then the result of conduc- 
tion and convection in both fluids, and is not constant 
above the onset of convection. The two-fluid interface is 
deformable, and its distance from the upper plate also 
changes the local heat flux. |l3|] In addition, buoyancy 
forces in the upper fluid may actively assist or impede 
convection in the lower fluid via surface stresses. |?],[l4],[l5] 
Recent experiments have found it necessary to treat the 
dynamics of both fluids in order the quantitatively ex- 
plain the data. fig] ] 

In order to gain insight into the more complex, but 
experimentally better posed, two-layer problem, we have 
undertaken an experimental study of convection in a sys- 
tem of two immiscible fluids sharing a deformable inter- 
face. To our knowledge, only one previous experiment 
has been done on this system, and it was restricted to 
locating the onset of convection, without visualization 
fuj . On the theoretical side, the two-layer problem has 
been the subject of a complete linear stability analysis 
frj, [l4|]T5| ] and some weakly nonlinear analyses (]]]24|] in 
certain limiting cases. The general problem, including 
surface deformations, should however be well within the 
range of modern weakly nonlinear theory and of numer- 
ical simulations. 



In our experiment, both fluids were liquids, and the 
conditions were such that the Marangoni forces were 
comparable to the buoyancy forces at onset. The di- 
mensionless heat flux across the interface is larger in our 
experiment than in previous studies. This is expected 
to have important effects in the nonlinear regime. |g2 24 
We used shadowgraph imaging to visualize the pattern 
and studied its secondary instabilities. We measured the 
mean wavenumber of the patterns, and compared them 
to the linear stability result, calculated for the full two- 
layer problem [p"5| . We could not resolve any hysteresis 
at the onset of convection. The first clear pattern we 
found above onset had a square planform. The mean 
wavenumber of this pattern is in good agreement with 
the critical value predicted from linear theory fl5|| . The 
square planform persisted up to e ~ 0.7, where it un- 
derwent a slightly hysteretic transition to a roll pattern. 
Here, e = (AT/AT C ) — 1, where AT is the temperature 
difference across both layers, and AT C is its value at the 
onset of convection. Up to the maximum value of e ~ 1.4 
we reached, the scale of the convection cells was deter- 
mined by the depth of the lower liquid. We compare our 
results with two recent nonlinear theories |^2],^4| and find 
qualitative agreement. 

The lower fluid was FC-75, a low-viscosity, perfluori- 
nated hydrocarbon fl6)| , while the upper was water. An 
important parameter is the depth fraction, L = d + /d, 
where d + (d) is the depth of the upper (lower) fluid. 
In this paper, we report results for L — 2.18 ± 0.04. 
The dimensionless heat flux across the interface, which 
is uniform below the onset of convection, is given by the 
Biot number B = A+/AL, where A + (A) is the thermal 
conductivity of the upper (lower) fluid 111 
to previous liquid/gas experiments HHjl 
which B ~ 0.1, in our experiment A~ > . 
of the total temperature drop AT falls across the thinner 
lower layer and 23 = 4.31±0.08at onset. 

Water, at a constant temperature T), bathed the top 
surface of a sapphire window which formed the top 
boundary of the cell. The base plate of the cell was 
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FIG. 1. Shadowgraph images of the patterns observed, (a) 
A disordered square planform at e = 0.48. (b) Near the transi- 
tion from squares to rolls at e = 0.78 , coexistence is observed, 
(c) Roll planform at e = 1.38. (d) A more ordered square pat- 
tern observed upon slowly decreasing e to 0.50, starting from 
roll pattern (c). 

maintained at a temperature T\ > Tq using a thin- 
film electric heater attached to its bottom surface. Both 
temperatures were feedback controlled to ±lmK. AT = 
T\ — Tq is the experimental control parameter. We es- 
timated the heat flux through the cell by monitoring 
the heater power. The sidewall of the cell was made of 
plastic and included a wedge-shaped radial fin to sta- 
bilize the two-fluid contact line. The round cell had 
thickness D = d + + d = 4.06 ± 0.01 mm, and radius 
r = 41.15 ±0.05 mm. The relevant aspect ratio is that of 
the lower fluid, T = r/d = 32.3. The uniformity of D was 
±8/im, checked interferometrically. The cell was leveled 
to within ±100/zrad, using an electronic bubble level. 

Both fluids wet all of the surfaces, so bridged config- 
urations of the interface were difficult to avoid. To pre- 
pare the two-layer configuration, the horizontal cell was 
completely filled with FC-75 which was then partially 
displaced by water slowly introduced above the fin. The 
interface was made level with the fin by visually elimi- 
nating the meniscus. The resulting metastable layering 
was somewhat imperfect due to the pinning of small de- 
formations of the contact line, and because the interface 
was difficult to see. The filling fraction L was determined 
destructively at the end of the experiment by tipping the 
cell on edge and measuring the interface position. 



The onset of convection and the nonlinear patterns 
were visualized by the shadowgraph method [pifl , using 
a beam reflected off the mirrored bottom plate of the 
cell. This beam passed twice through each fluid and twice 
through the interface. The shadowgraph signal therefore 
contained contrast due to temperature gradients in both 
fluids, plus a component due to deflections at the de- 
formed interface. The latter effect is likely to be small, 
since the two fluids are nearly index matched. We ob- 
served that the scale of the pattern in the upper layer was 
slaved to that of the lower. Thus, it was not necessary to 
disentangle the various contributions to the shadowgraph 
images in order to understand the basic planforms. 

A typical run consisted of slowly increasing and de- 
creasing AT, while recording images and measuring the 
heat flux. The onset of convection occurred at AT C = 
0.999 ± 0.025°C. Convection entered from one side of 
cell as a highly disordered pattern, and eventually filled 
the cell. We did not observe hysteresis in the onset, pre- 
sumably because of the rounding effect due to nonunifor- 
mities in the cell. The first clear patterns that emerged 
had a square planform. Square patterns have been ob- 
served previously in Marangoni convection Jll| , but only 
at much higher e as a secondary bifurcation from a hexag- 
onal pattern. Here, our observations are consistent with a 
square planform emerging at onset, but, due to rounding 
effects, we cannot exclude the possibility that a hexago- 
nal pattern exists in a narrow range close to onset. 

Typical patterns are shown in Figure 0. For increasing 
e, a patchy, time-dependent square planform, as shown 
in Figure 0(a) was found between onset and e sr — 0.70 ± 
0.01 where a slightly hysteretic transition to rolls was 
observed. This transition could easily be located by a 
significant increase in the heat flux through the cell. Near 
the transition, we find a dynamical coexistence of squares 
and rolls, as shown in Figure 0(b). The rolls appeared 
where squares merged along shared edges, and visa versa. 
At values of e > e sr , the pattern was dominated by rolls, 
shown in Figure 0(c), with squares appearing only at 
grain boundaries. The more ordered roll patterns were 
very slowly time dependent. The finned sidewalls had 
a rather weak orienting effect on the rolls. As e was 
decreased, an ordered square pattern re-emerged at e rs — 
0.67 ± 0.02 from the pre-existing rolls. The hysteresis 
e sr — e rs = 0.03 was just barely resolvable from the heat 
flux data. Figure 0(d) shows such a square planform. In 
all cases, the patterns had a scale ~ d, the dimension of 
the thinner lower layer. 

The characteristic timescale for convection in the lower 
fluid is the vertical thermal diffusion time t v — d /k — 
47s, where n is the thermal diffusivity. We ramped the 
temperature slowly compared to t v , but it was impracti- 
cal to increase it more slowly than the much-longer hor- 
izontal diffusion time = T 2 t v — 14hrs. Thus, some 
of the disorder we observe on increasing e may be due 
to the finite ramp rate. However, we also performed ex- 
periments in which e was ramped very quickly to a value 
< e sr , and then held constant for ~ 6^. The resulting 
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k c = 2.38 



FIG. 2. The mean wavenumber (fc), showing the linear sta- 
bility boundary from Ref. |l5| (solid line). Bars indicate the 
spread in (k) due to disorder in the pattern. We find a square 
planform (square symbols) for e < 0.7, and rolls (circular 
symbols) otherwise. Open (solid) symbols indicate increasing 
(decreasing) e. 



square pattern did not anneal significantly and re- 
mained time-dependent. Thus, the disorder may be par- 
tially dynamical in origin. 

We determined the wavenumber of the patterns from 
Fourier analysis. We found the azimuthally averaged 
structure function of a 256 2 pixel region of the centre 
of the cell. The peak in the structure function gave the 
mean wavenumber (k). Fig. || compares (k) with the 
results of linear stability theory |l5f| . The dimension- 
less wavenumber k is scaled by the depth of the bottom 
fluid d. The value of (k) at small e agrees with the crit- 
ical value k c — 2.38 from the linear theory. The more 
ordered square patterns observed for decreasing e have 
slightly larger (k) than those of increasing e. 

To evaluate the linear stability boundary shown in 
Fig. |^, all of the material parameters are needed. These 
are, for the lower fluid, the kinematic viscosity u, the 
density p, and k = A/pCp, where Cp is the heat ca- 
pacity, as well as the corresponding quantities for the 
upper fluid and the interfacial surface tension a. All 
of these are assumed to be independent of T, accord- 
ing to the usual Boussinesq approximation, except a and 
the two densities, which are linearly temperature depen- 
dent. Thus, we also require the thermal expansion coeffi- 



cient P = (1/ p)dp/dT for each fluid and the temperature 
derivative of the interfacial surface tension 7 = —da/dT. 
The negative sign is included so that 7 > 0. All of the 
parameters arc known, except er and 7. Using Antonow's 
rule , we estimated a to be equal to the difference of 
the surface tensions of the two fluids, measured against 
air. This approximation suffices, because the neutral sta- 
bility boundary is very insensitive to a, which only en- 
ters into the small surface deformations JXSj ] . On the 
other hand, the remaining parameter 7 is crucial to the 
Marangoni mechanism. We fixed 7 by requiring that the 
linear theory reproduce the correct measured AT C . The 
result is 7 = 0.047±0.003 dynes/cm K. This is consistent 
to within a factor of two with the difference of the known 
values of 7 for each fluid, measured against air (i.e. with 
Antonow's rule jl8| applied to da/dT). 

One- and two-layer theories of Benard-Marangoni con- 
vection traditionally scale the problem using the temper- 
ature difference AT L across the lower fluid. Below onset, 

AT L = (1 + AL/A+)- 1 AT = (1 + 1/B) _1 AT. (0.1) 

This scaling is only approximate above onset, where con- 
vection effectively increases the average thermal conduc- 
tivity of each layer, and also introduces temperature vari- 
ations on the interface 0,0- The often- used assump- 
tion of a constant B is problematical above onset for the 
same reasons. [ [l7||2"cfl However, in order to make contact 
with theo ries using this scaling, it is customary to adopt 
Eqn. 0.1 even above onset. M The seriousness of this 



approximation can be appreciated from heat flux mea- 
surements. The total effective thermal conductivity 
of the cell at the highest e was only about 8% larger 
than that below ons et. This sets an upper bound on the 



validity of Eqn. 3.1 



Assuming Eqn. 0.1, one can define the Rayleigh and 



Marangoni numbers 7Z and yVf of the lower layer as fol- 
lows; 
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where g is the acceleration due to gravity. TZ and M. are 
not independent, so it is convenient to define two related 
parameters |23j], 
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Here, TZq (A^o) is the critical value of 1Z (A4) in the ab- 
sence of Marangoni (buoyancy) forces. In general, in the 
two-layer problem TZq and A^o will differ from the usual 
critical values for any one-layer model of the lower fluid. 
They can be determined from the full linear theory of 
the coupled, two-layer system, as described below, a is 
fixed by the choice of fluids and the depth d, and can 
be considered as the fraction of the total forcing due to 
buoyancy, while A is proportional to the total tempera- 
ture difference AT. 
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FIG. 3. The Marangoni number M vs. the Rayleigh num- 
ber 1Z for the lower fluid, showing the neutral stability bound- 
ary (thick solid line) from the complete two-layer analysis of 
Ref. [jl5| . Open (solid) symbols show data for increasing (de- 
creasing) e. Square (circular) symbols represent points where 
the square (roll) planform was observed. 

Fig. U shows a plot of M. vs. 1Z. The linear stability 
boundary for the coupled, two-layer system at constant L 
separates the conduction regime near the origin from the 
convecting regime. Its intercepts are the critical values 
TZo and Mo of each number in the absence of the comple- 
mentary effect. The boundary is calculated by holding 
L and the material parameters constant and varying d. 
One point on this line corresponds to the experimental 
value of d. a varies from on the M axis to 1 on the 
7Z axis. Above the line a = 1/2, Marangoni forces domi- 
nate; below, buoyancy. For our experimental conditions, 
a = 0.58, so tha t the two effects are nearly equal. As- 
suming Eqn. 0.1 holds, increasing A oc e corresponds to 
moving outward along the radial line a = 0.58. The 
symbols show the values of M and 1Z that we reached, 
and the point where the planform changed from squares 
to rolls. It would be interesting to map the stability 
boundaries of various patterns in the (1Z,M) plane by 
experimentally varying a and A independently. This can 
be done in principle by varying e and d, keeping L fixed. 

Two recent nonlinear theories are relevant to our obser- 
vations ^H,^4) , although neither treats the full two-layer 
problem. Regnier et. al. included both surface ten- 
sion and buoyancy in a one-layer model with no surface 
deformations and a constant B boundary condition on 



the upper surface. They derived a set of amplitude equa- 
tions which show various regimes of hexagonal, square 
and roll planforms above onset. Their main result is that 
there exists a critical value a c above which a hexagonal 
pattern at onset gives way to squares for a < a c or to 
rolls for a > a c . Unfortunately, they only present results 
for B < 2, where our B ~ 4. The trends however suggest 
that as B is increased in their model, a c — > 1, favouring 
squares. In the same limit, the range above onset over 
which the hexagonal pattern exists becomes very narrow. 
These two predictions are in qualitative agreement with 
our observations at low e. 

Golovin et. al. |24| treated the full two-layer prob- 
lem, but included Marangoni forces only. They derived 
amplitude equations, including the effect of surface defor- 
mations, but neglecting stresses on the interface on the 
assumption that the upper fluid was a gas. Interestingly, 
when the heat transfer through the upper fluid (i.e. B) 
is increased, a square pattern is predicted at onset. This 
is also consistent with our observations. Neither of these 
nonlinear theories predicts the square to roll transition 
we observed, however. Nor do the amplitude equations 
contain the gradient terms that would be necessary to 
describe disordered patterns of the kind we observe. 

In summary, we have performed experiments on ther- 
mal convection in a system of two layered liquids under 
conditions where both buoyancy and Marangoni forces 
are comparable. We found a square planform of convec- 
tion just above onset, which underwent a transition to a 
roll pattern at higher control parameter. We compared 
the wavenumber of the patterns to the linear stability 
boundary for the full two-layer coupled problem. We 
found quantitative agreement between the wavenumber 
of the pattern just above onset and the critical wavenum- 
ber k c from linear theory. Two nonlinear theories, orig- 
inally proposed to explain gas over liquid experiments, 
can be applied to our case. They qualitatively predict 
the square patterns found near onset, but not the ob- 
served transition to rolls. 
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